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A previous paper derived the solution for the drag-free motion of an artificial satellite 
in the gravitational field of an oblate planet. The corresponding potential, expressed in 
oblate spheroidal coordinates, leads to separability and represents the even zonal harmonics 
exactly through the second, for any oblate planet, and approximately through the fourth, 
in the case of the earth. 

The previous paper contained a restriction on the orbital inclination 7, viz, /c<^<Cl80° 
— Ic, where Ic might be as large as 1°54' for an orbit sufficiently close to the earth. The 
present paper removes this restriction and shows that many of the formulae for the periodic 
terms may be simplified, when the orbit is equatorial or almost so. The results agree with 
those obtained by a direct two-dimensional solution, when the orbit is purely equatorial. 

1. Introduction 

This paper is a sequel to a recent paper,^ concerning an accurate intermediary orbit for 
satellite astronomy, and will accordingly follow the notation thereof. It there followed that if 

X^Ji/62<l, (1) 

all the p-integrals are expressible in terms of rapidly converging series involving products of 
Legendre polynomials with arguments X and (l—e^)~'K Condition (1) is equivalent to a 
restriction on the orbital inclination /, viz, 

/,<7<180°-/,, (2) 

where, to the first order in k, 

tun' I,=k^{rjpyj2. (3) 

For the earth e72 = 0. 00108, so that for orbits so close that ])^rej 

/c=l°54'. (4) 

I imposed the condition X<^ 1 in order that 

(l+^p"-i+5p-2)-^ = (l-2X/i+A2)-^ (5) 

should be a generating function for the Legendre polynomials P,XX). It now appears that 
such a restriction is unnecessary. Thus 

{l-2\h + h'y'^ = ^h''Pn{\) (6) 

even when X^ 1, provided only that ^ 

h<\-(\'-l)h (7) 

We need show only that (7) is always satisfied, in order that all the results of (A) hold for all 
orbital inclinations. The only changes will be a few simplifications, especially for the cases 
7=0° or 180° of purely equatorial orbits. 



1 This work was supported by the U.S. Air Force, through the Office of Scientific Kesearch of the Air Research and Development Command. 

2 J. P. Vinti, J. Research NBS, 65B (Math, and Math. Phys.), 169-201 (1961), hereafter referred to as (A). Any reference in the present paper 
to an equation with a decimal number, such as (5.30), denotes an equation in (A). 

3 E. W. Hobson, The Theory of Spherical and Ellipsoidal Harmonics, p. 15, Cambridge University Press, Cambridge 1931. 



2. The p-Integrals 

To show that (7) is ahva3"s satisfied Avhen X^ 1, note that 

bi hi 62 



-A -^^' (8) 

P 02 p 

By (4.12), (4.13), and (5.6) it follows that 61 and 62 are both real and non-negative for c^<^(ip, 
a relation that always holds for satellite orbits. Thus (7) is satisfied if and onl3" if 

&i/p<^(X), (9) 

where 

g(X)^X'-\{\'-l)K (10) 

We shall next show that 

1/2<^(X)^1 (X^l). (11) 

To show that ^(X;^l, note that for X^l, we have 1>1 — X'^^O, so that {l — \~^)'^l 
-X-2, i-(i_x-2)i^x-^ and finally X2-X(X2-1)^^1. Thus ^(X) ^1. 
To show that (/(X)>l/2, note that for X^ 1 

0^1-X-2<(l-JX-2)2, (12) 

so that 

(i-x-2)Ki-ix-2 (13) 

or 

l_(l_X-^)i>iX-l (14) 

Then 

(7(X)=V-X(X2-l)i>l/2, (15) 

as was to be shown. 

From (9) and (11) it follows that when X^l the condition 6i/p<Cl is necessary for the 
validity of (7) and that the condition 6i/p<l/2 is sufficient for its validity. 
To show that 6i/p<l/2 for all bound * orbits, note that p^pi, so that 

6i/p^6,/p,. (16) 

From the relation c^=kp' and the relation pi=a(l-e),e^l, for a bound orbit, we then find 
from (3.25), (4.12), and (5.6) that 



b^__ Hl+e)il-vl)ll-kvlil-e')] 



(17) 



(18) 



p, [l-k{l-e')][l-kr,l{i-e')]+4kr,l 
Here the numerator ^^(1+e) and the denominator ^(l—k)^, so that 

h^ kjl+e) ^ 2k 

pr {i-kr= {l-kf 

a function monotonic in k for 0<Z:<1 and <l/2 whenever ^<3 — 2y^(==0.l72). 
For the earth k<C0M108, so that 

&i/pi<0.00216<l/2. (19) 

From (16) and (19) it follows that for the earth 

^>i/p<l/2. (20) 



* Past participle of tho verb "to bind," taken from the terminology of atomic theory. 



This is the condition sufficient for the vahdit}" of (7) and thus of the expansion (6), for all 

All the developments of (A) then liold for the p- integrals, in particular (5.30) through 
(5.33) for the integral i?i, (5.35) through (5.40) for B2, and (5.60) through (5.65) for E^, where 
the Dm!^ are again given by (5.50) and (5.53). 

To show the rapid convergence of the various series that there occur, w^e note first that since 

Pn{^)=T^-' V [\+{\'-iy'\^OSXfdx (21) 

and since X = 6i/62^1, we have 

{hlvYPniA-i^-' J^^ [h,r'+P~\b\-h'^"' cos xYdx, (22) 



Then 



|(^)>.(X)|^g+^ (1-X-T-J^(f J- (23) 



From (5.14) and (5.34) the series Si and S2 that occur in the expressions for the p-integrals 
Ri and i?o are 

S,^f^CiyPn{\) f (1 + 6 cos^)"-^.^^, (i=l,2) (24) 

ii=iiiW / Jo 

where % = 2 and 712= 0. Thus 

\Sj\^v±(^-^y{l + er-j (25) 

..S(|f^i-^^)''.<t)Mj[fa.^)T 

^-^^ ■ (i=l,2). (27) 

l-|l(l+.) 

But p=a(l—e^) = pi(l + e), so that by (18) 

b,p-'Sk{l-k)-' (28) 

and 

2bi(l+e)p-' S2k(l+e)(l-k)-'^4.k{l-k)~', (29) 

where 4^(1— ^)~^<0. 0043 for the earth. It follows that the series for 7?i and 7?2, and thus the 
series for the secular coefficients Ai and A2, converge absolutely and more rapidly than a 
geometric series of common ratio 0.0043. 

By (5.49), (5.50), and (5.53) the series 5^3 that occurs in the expression for the p-integral 
i?3 is 

S,=Se+So, (30) 

Se^S A. (\l + e COS xy-+'dx (31) 

71 = Jo 

So^J:D,„+, Cn + e cos j^)2"+3 dx (32) 

71 = Jo 



A.=S(-l)"-^(c/^)'"-''(&2/^)'^^2;(X) (33) 

j=0 



From (33) and (23) and the relation c^=kp^, it follows that 

j=0 j=0 



= i-(2bjcy 

and from (34) and (23) that 

\D2n+,\^k^{2bJp) ± {2h,lcy^^k\2b,lp) ± (2h,lcr 
•■=0 ;=o 



Thus, by (31) and (35) 



and by (32) and (36) that 



'^^l = l-(26i/c)^^o^ ^ ^ ' 

(1 + e)^^ 
= [l-(26,/c)^][l-i(l+e)^] 

(l + e)%(26,/y) 
= [l-(26:/c)^][l-A:(l + e)2]' 



(35) 



<i!(2V|?L. (36) 



(37) 



(38) 



Then, by (30), (37), and (38), we have 

I l^ (l + e)Ml + (l + ^)(26,/y)] .„„. 

\^'\= [i-{2b^icy\[i-k{i+ey\ ^-^^^ 

But p=pi(l+e), so that by (18) 

2b,/p^2kil-k)-' (40) 

and 

2b,/c^i2b,/p)k-"'S2k"Kl-k)-'. (41) 

Thus 

|o|^ (l + g)Ml+2A:(l+e)(l-Ar)-^ ] .,„, 

'*^''= [l-4^(l-^)-«][l-t(l + e)^]' ^* ^ 

The series for Sa thus converges absolutelj^ and more rapidly than the power series in ^ of the 
function on the right side of (42), where ^<^0.0011 for the earth. On replacing v by x, we can 
then say the same thing about the secular coefficient A3. 

3. Simplification of the p-Coefficients When X^ 1 

From X= 61/62 it follows that 

bi=bi\-'^bi (X^l). (43) 

Since bi = 0{k), it then follows that 62 is also of order k when X^ 1. This fact enables us to sim- 
plify the coefficients Am and Aq,,, which are needed only to 0{k^), and the coefficients A^n, 
which are needed only to 0(k). 
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Thus (5.32) and (5.33) lead to 

An=0{l^) A,,= 0{¥), (44) 

(5.39) and (5.40) lead to 

^23= 0(F) ^24=0(P), (45) 

and (5.37) and (5.38) lead to 

^2i=(l-g2)ip-ie[6,2>-i + (3-X-^)Fcos*/] (46) 

A22={l-e')^p-' ~ (3-X-^)F cos*/. (47) 

Here we have used 

b,p-'=k cos' I +Oik') (48) 

in the terms involving^&j p~'. 

Similarly, for the coefScients A^n, which are needed only to 0(k)j we obtain from (5.62) 
through (5.65) 

A,,= {l-e')^p-'e r2+(3+| e'") k cos^/- (4+36^)^1 (49) 

A,,= {\-e^)h^-H' [i+^ k cos-^ /~(^+f) ^] (50) 

A,,= {l-eynr'e' [^ cos^^-^] (51) 

Au==-{\-e'y^p-H'kl'^2. (52) 

For X^ 1 we need also to rewrite (5.31), (5.36), (5.50), and (5.53) for tlie secular coefficients 
Ai, A2, D2n, and jDsw+i, which all contam terms of the form (b2/2^)^PmW - As sin / approaches 
zero, 62 also approaches zero and X becomes mfinite, so that such a term takes the indeterminate 
form zero times infinity. To remove tliis indeterminacy, note tliat 



so that 



*^^^^*.^^X-', (53) 

(|)"'7^„(X)=(|y'x-'«P„,(X) (54) 

=(^)'"7.',„(X-). (55) 

Here Rm{x) is a function that has already appeared in (A), viz, 

R,,{x)^X-^Prr,(x") (56) 

a polynomial of degree [m/2] in x'. 

To determine (b2/p)'^PmW for sin 1=0, first write 



n.(A)-^ 2-i!(m-i)!(m-2i)!' 



(57) 



so that 

i?;.(X-^)^X— P^(X) = 1 Cm-0, 1) 

__(2m)J_ M (-l)^(2m-2j)!X-^^' 
~T'(m\y'^pi 2"^ j\(m-j)l(m-2j)\ 
Thus 



(m=2, 3, 4, . . .) (58) 



^U0)=2^^; (^-0,1,2,3, . . .) (59) 



620565—62 2 



so that by (55) and (59) 

(y^)-PJX)=|^;(^)", (sin/=0) . (60) 

where bio is the value of bi for T/o^^siii 1=0. 

From (5.31), (5.36), (55), and (60) it then follows that 

Ai={l-er^p f: {b,lpYRn{\-')Rn-2Ul^') (X^ 1) (61) 

71 = 2 

= (l-e^)¥ ± ||] (^y iZ„-2(VI^) (sin 7=0) (62) 

A,= {l-e'yr'T.{hlpYRn{\-')Rn(d^^^'-) (X^l) (63) 

= (l-.^)i^-> g Hi (^)"i?„(v'l^) (sin 7=0) (64) 

Note that A^=0{¥) and A=0(^'>). 
Then, from (5.61) 

yl3=(l-0¥-'f: ^««»+2U/l^^), (65) 



w=0 



where, by (5.50), (5.53), (55), (60), and the relation c^=kp^, we find 

j=0 



=2] i-ir-'k-^{bJpy^R2j{\-') (X^l) (66) 



;=0 



7=0 

=Z: (-l)"-^A:"-^(^/p)«+'iZ2,+i(X-0 (X^l) (68) 

4. Simplification of Other Coefficients When X^ 1 

To simplify the 77-coefficients when X^ 1, we must first note that from page (176) of (A) 
we have 

bi^kp cos^ / 62 -Jc^p sin /, (70) 

so that 

X = &i/62«^*cos2/csc/ (71) 

or 

r;o=sin I^kh-' cos^ /. (72) 

Thus 

vo=0(ki) (x^i). (73) 

10 



Also, since r? 2=0{k^), hj (3.42), it rollovvs that 

q^Vo/v2=0(k) (X^l). (74) 

There also arises an mdeterminacy in the 77-coefficients, from the quantity (al—af)~^rioy 
which takes the form infinity times zero as 77o = sin / approaches zero. To remove tliis in- 
determinacy, use (4.15a), viz, 

as=a2(l--^Y' cos /, (75) 

to find 

(ai-ai)-"^'„o=«.-' fl +-4- COS^ /)"'"• (76) 

Then, from (8.27) and (7.23) 

4^,=27rvAt+Pi+l32a^Ka+b,+AM2'] (77) 

2Trv2=a2(l+—cos'i)' A^B^' (a+fti+^i+c'r/^i^Bii^^"')"'. (78) 

From (8.33) 

^0= {-2a,)-'''a, (l+— cos^ /Y''yl,5,-^z;o (79) 

From (8.37) 
M,= {a+h,)-' [-{A,^-c'rilA,BJ3^') v,+~ {-2a;)^a^'(l^-~vo^ lY"^ y{i^ (2^.+2v^o)], (80) 

of order P. Then, by (8.39) 

E,= [l-e' cos {M,+E,)]-'M,, (81) 

since M\=0{¥), Then E, and v, are botli of order k\ Also, by (8.40), 
^1- {-2a,) -^ a, (l +— cos^/Y B^'[A2V,+A2, sin (7l/,+z;o) 

+^22sin (2M,+2^'o)]+|52-^shi (2^,+2^o). (82) 

By (8.45), M2 and thus E2 and ?;2 <ii'0 of order B and, by (8.48), 1/^2 is also of order P. Thus 
all the second-order periodic terms of (A) becojue of tlie third order and thus neghgible, when 
X^l. 

B}^ (8.50) the right ascension ip becomes 

^=/33+a3«2~^(l+-^COs2/y'[(l-r7g)-Ul-r72-')-^X + 53^] 

-c'a, (-2aO-^(^3?^+g ^.sin/^z;). (83) 

Here we have dropped the periodic term (3/32) 'nl'n2^ sin 2\l/ of (8.50), since it is of order P for 
X^l. 

5. Summary for X^ 1 

For an ahnost equatorial orbit, corresponchng to X^l, tlie right ascension ^ is given by 
(83). The st)heroidal coordinates p and 77 are given by 

p=a{l — e cos£) = (l+e cos 'v)~^p 

r, = rjo sin xly. 
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Here the expressions 

E=Ms+Eo+E, v=Ms+Vo+v, ^-\^.+fo+\^i 

are sufficiently accurate to give the secular terms exactly and the periodic terms correctly 
through order k\ provided that Ms is calculated by (8.24), ^, by (77), Eo by (8.31), Ei by (80) 
and (81), Vq and Vi by the anomaly relations (8.1), xpo by (79), and ^i by (82). 

6. The Case of a Purely Equatorial Orbit, 7=0° or 180° 

For 7=0° or 180° we have 770=0, cos' 1=1, x^^ by (6.51), y„ = ^ by (6.66), 

53-l-(l-r72-')-^ (84) 

by (6.65), and |a3l = a2, so that 

«3/«2=sgn 0:3. (85) 

Then, by (83), 

^=/?3+(sgna3) (l+-^\" yP-c'a,{-2a,)-'^ (^3^+ 2 A« sin m^\ (86) 

Also, by (8.24) with 770=0, 

M.= {^-2a)Ka+h,+A,)-\t+^,) (87) 

and by (77) and (78), with 770=0, cos' ^=1, ^^^i^^l ^2===1, 

yP,=pJl+^\^+a,AJl+^\\a+^^^ (88) 

Thus, by (87) and (88), 

^^=^'(i+;r^y+^'(i+;r^y (-2^.0-* a^m.. (89) 

\ <^oPo/ \ (^vPo/ 

Then, by (79) and (82), with 770=0, 

^o=^'(i+;rJry (-2«i)"*^2^'o (90) 

,^i=a,(l+-^Y (-2aO-*^^2^'l+S^2nsin(niK+w■o)l• (91) 

\ CioPo/ L. n=l J 

Addition of (89) through (91) then gives 

^^^' (l +;rJry +^' (^ +^1 ^-2^^) ~* l^^v+i: A^, sm(nMs+nvo)']' (92) 

V (loi^o/ \ (^oPo/ L n=l J 

Since ?;=Ms+z?o+^b where ?;i is of order A:^, and since A21 and A22 are of orders k and /j^ respec- 
tively, it follows that 

2 2 

2 ^2^ sin(7i,M,+m^o) = X1 ^2« sin nv (93) 

to order k^. Thus, to order k^, 

On inserting (94) into (86), we then find 

(t>^p3 f & sgn a3+Q:3(— 2ai)"* (A2V+^ ^2. sin nvj—c''as{-2a;)-i (a^v+^ A^n sin nv\ (95) 
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It is a simple exercise to check these results for a purely equatorial orbit. To do so, let 
X and Y be the usual Cartesian coordinates, define p and <i> by 

X+ir=(p2+c2)* exp ^>, (96) 

write down the kinetic energy |(X^+F^) in terms of p and </> and their time derivatives, write 
the potential as — mp~S construct the Hamiltonian, and then write down the Hamilton-Jacobi 
equation. Separate the latter, to obtain the solution 

W^a^± f {p'+cY' F{p)i dp, (97) 

where 

F(p)^c'al + (p'+c')(-al+2np+2a,p'), (98) 

The kinetic equations are then 

bW Cp 1 

^=t+l3,= ± p'F{p)-Ulp (99) 

|^=03=<^T«3 f (p'+cT' p'F{p)-^ dp (100) 

=<^Ta3 ( ' Hp)-^ dp±c'a-, i " (p2+C^)-' F{p)-i dp. (101) 

J Pi J Pi 

On then following the procedure in (A), we find that p is given by just the results of section 5 
of the present paper, with 770 placed equal to ya'vo, and that is givcMi by 

0=^3+«3(-2a:i)-^ (AoV+^^ A2n sin nv)-c'a, (-2^0"^ (a^'+^^ An sin m\ (102) 

In order to reconcile (102) and (95), observe that for an equatorial orbit ^2 ^i^d fe ^re not 
defined separately, but enter in the combination 03+132 sgn 0:3. Accordingly, the definition 

0z=0z+02 Sgn ^3 (103) 

will bring (102) and (95) into perfect agrcH^nuMit. To understand this definition better, we 
reason as follows. 

In the equatorial plane /'^=p^+c^, so that r is at mininiuin whenever p is at minimum. 
That is, the satellite is at perigee wlienever v=2TrT, r=0, 1,2,.... By (102) the right ascension 
</)p thus changes value from one perigee to the next. For an equatorial orbit about a planet of 
zero oblateness, however, the coefficient of v in (102), viz, as{ — 2ai)~^A2—c^As)j would reduce 
to sgn a^, so that in such a limiting case 

c/>,=/5^+2^r sgn ^3,(^-0, 1, 2, . . .) (104) 

so that the actual position of perigee would remain fixed and (3^ would be its right ascension. 
But for the case of a nonequatorial orbit around a planet of zero oblateness 0s is simply 
the right ascension of the ascending node and 02 is the argument co of perigee. Thus for the 
limiting case of a purely equatorial orbit about a planet of zero oblateness (103) would take 
the expected form 

R.A. of perigee=l^±co, (105) 

where the sign would be plus for a direct orbit and minus for a retrograde orbit. Thus the 
definition (103) is reasonable. 

(Paper 66BI-68) 
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